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ABSTRACT 

Context. Most PG 1159 stars are hydrogen deficient post-AGB stars on their way to the white dwarf cooling sequence. However, a 
fraction of them show small amounts of atmospheric hydrogen and are known as the hybrid PG 1 159 stars. Both the normal and hybrid 
PG 1 159 stars may reside at the centre of a Planetary Nebula. Some PG 1 159 stars exhibit non-radial ^-mode pulsations, providing us 
with a means to constrain their fundamental physical parameters, internal structure, and evolutionary status; asteroseismology is the 
technique of using global stellar pulsations to probe the star's interior. 

Aims. Model calculations demonstrate that the PG 1159 instability strip is driven by the partial ionization of carbon and oxygen 
through the ^c-mechanism. The range of excited frequencies in models is sensitive to the details of the chemical composition in 
the driving zone. Hydrogen content plays a crucial role in deciding the range of excited frequencies for models of hybrid PG 1159 
pulsators. The range of periods observed in all PG 1159 pulsators can then be compared against the range predicted by theoretical 
models, which are based on observations of their chemical composition and evolutionary status. Our aim is to test the validity of the 
predicted range of pulsation periods for the hybrid PG 1 159 variables, Abell43 and NGC 7094. 

Methods. To achieve this goal, we acquired continuous high speed CCD photometry of the hybrid PG 1159 star Abell43 using the 
2.6 m Nordic Optical Telescope and the 3.5 m telescope at Apache Point Observatory. We also observed NGC 7094, a hybrid PG 1 159 
star with atmospheric parameters similar to Abell43. We analyzed our data to obtain the pulsation spectra of the two stars. 
Results. We detected six significant pulsations in the light curves of Abell43 with periods between 2380 s and 6075 s. We were also 
fortunate to detect low amplitude pulsations in our NGC 7094 data with periods between 2000 s and 5000 s. The observed range of 
periods for both stars is consistent with the theoretical predictions. 

Conclusions. The range of periods observed in both hybrid PG 1 159 stars, Abell 43 and NGC 7094, agree with the theoretical model 
that the pulsations are driven by the /c-mechanism induced by the partial ionization of carbon and oxygen. 

Key words, stars: planetary system - stars: oscillations - stars: individual: Abell 43, NGC 7094 central star. 

1. Introduction Low carbon abundance determinations by Dreizler et al. 

(119951) had initially led Quirion et al. (120041) to rule out pul- 

Hydrogen deficient hot white dwarf stars constitute the PG 1 159 sations in Abell 43 based on their model calculations. But the 

class, and are often found at the centre of planetary nebu- recent high carbon abundance measurements led Quirion et al. 

lae. Abell 43 (=WD 1751 + 10) is one of the four known hybrid ( 120051 ) to re-evaluate their previous conclusions about Abell 43. 

PG 1159 stars, which show evidence of a non-negligible hydro- They found that the revised models suggested that Abell 43 

gen abundance in their atmospheres (Napiwotzki & Schonberner should exhibit { - 1 g-mode pulsations, driven by the classical 

11991) . Miksa et al. (,2002j recently pubUshed the following at- ^--mechanism due to the ionization of carbon K-she ll electrons, 

mospheric parameters for Abell 43: Teff = llOkK, log g = 5.1, in the range of ~2604-5529s. Corsico et al. (120061 ) carried out 

X(H) = 0.35, X(He) = 0.42 , and X(C) = 0.23 in mass fractions, nonadiabatic pulsation calculations for a grid of evolutionary 

Werner & Herwig ( 120061 ) found an even higher carbon abun- PG 1159 models. They predicted that Abell 43 should exhibit pe- 

dance from HST ultra-violet observations of Abell 43. riods in the range of -2600-7000 s adopting a O.53OM0 stellar 

mass evolutionary track. 

Send offprint requests to: J.-E. Solheim CiarduUo & Bond (11996D were the first to report possi- 

* Based on observations made with the Nordic Optical Telescope, op- ^le variations in the light curve of Abell 43 with a period of 

erated on the island of La Palma jointly by Denmark, Finland, Iceland, about 2473 s. Schuh et al. (|2000|) detected a longer period of 

Norway, and Sweden, in the Spanish Observatorio del Roque de los 5500 s. They also suggested that the -2500 s period reported 

Muchachos of the Institute de Astrofisica de Canarias. by Ciardullo & Bond ( 1996 ) could simply be an artifact of the 

** Hubble fellow 5500 s period caused by the short duration of the observations. 
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Fig. 1. Light curves of Abell 43 shown for the three long NOT 
runs. The solid lines are synthetic light curves constructed from 
the parameters listed in Table 2, and shifted downward by 0.005 
units for clarity. 

Table 1. Log of observations of Abell 43 



Date 


Start 


Length 


Cycle 


Data 


Observer 




(UT) 


(hrs) 


time(s) 


points 




2005-06-29 


21:28 


4.40 


30 


541 


JES" 


2005-07-01 


03:43 


1.56 


20 


283 


JES" 


2005-07-01 


09:14 


1.54 


37 


51 


ASM^ 


2005-07-01 


21:08 


7.71 


20 


1399 


JES~ 


2005-07-02 


21:07 


7.32 


20 


1319 


JES 


2005-07-10 


07:24 


1.71 


41 


110 


ASM 



J. - E. Solheim 
A. S. Mukadam 



Vauclair et al. (120051) discovered multiple modes in the pulsation 
spectrum of Abell 43, and reported periods at 2600 s and 3035 s. 
Their observations were based on three individual runs, amount- 
ing to a net observation time of 7.5 hr at the 2.6 m Nordic Optical 
Telescope (NOT). 

In this paper, we report the results from 24 hours of time- 
series photometry acquired on Abell 43 in 2005 with the 2.6 m 
NOT and the 3.5 m telescope at Apache Point Observatory 
(APO). We also report pulsation periods detected in the hy- 
brid PG1159 star NGC7094 (= WD 2134-H125), which is al- 
most Uke a spectroscopic twin of AbeU43. 

2. Observations 

We observed Abell 43 for four nights, from the 29th of June to 
the 2nd of July, 2005 using the 2.6 m NOT in addition to two 
short runs with the 3.5 m telescope at APO on the 1st and 10th 
of July, 2005. The log of observations is given in Table [T] 

We conducted observations at the NOT using the Andalucia 
Faint Object Spectrograph and Camera (ALFOSC), which is 
equipped with a thinned 2048x2048 E2V CCD 42-40 chip. 
We observed with the broadband blue W92 filter, which has a 
FWHM of 275 nm centered at 550 nm; this enables us to gather 
a fair amount of flux from the target and minimize the sky con- 
tribution at the same time. We acquired suitable flat fields each 
night just after sunset and determined the dark current for each 
individual image using the overscan region. We controlled the 
ALFOSC data acquisition using the software interface tcpcom. 




Fig. 2. Pulsation spectrum of the light curve of the NOT run from 
the 2nd of July, 2005. The insert is a blow up of the low fre- 
quency region, that shows the evident lack of excited frequencies 
above 0.5 mHz. 



defining multiple windows to achieve fast readouts (0stensen 
i2000)- We achieved a net readout time of 4.8 to 5.8 s for the 
six separate windows around the target star, three comparison 
stars, and two blank fields for determining the simultaneous 
sky brightness. We chose a time resolution of 30 s during our 
first night at the NOT telescope, and our exposure time became 
23.9 s. We reduced the cycle time to 20 s on subsequent nights, 
and achieved exposure times of about 14.1, 14.8, and 14.9 s for 
slightly diff'erent windowing patterns. With real time process- 
ing, we were able to display light curves of the raw and sky- 
subtracted data during acquisition. 

We followed the standard basic reduction procedure, sub- 
tracting a dark frame from each image and then flat-fielding it. 
After these preliminary steps, we computed light curves using 
aperture photometry for several aperture sizes. We selected the 
best light curve based on the highest S/N ratio. For the observa- 
tions outlined above, we obtained the best result for an aperture 
diameter of 24 pixels corresponding to 4.6 arcseconds on the sky. 
We corrected for transparency variations and differential extinc- 
tion by dividing the target star light curve with the average light 
curve of the comparison stars. In Fig. [T] we show the light curves 
from the three longest NOT runs. 

The APO observations were acquired with the Seaver 
Prototype Imaging camera (SPIcam) and a Johnson B filter We 
used 3x3 binning and windowing in order to reduce the read 
out time to about 25-35 s. The instrument has a read noise of 
5 electrons RMS, a gain of 1 e/ADU, and a plate scale of 0.42 
arcsec/pixel with 3x3 binning. We used a standard IRAF reduc- 
tion to extract sky-subtracted light curves from the CCD frames 
using weighted circular aperture photometry (O'Donoghue et al. 
120001) . 

3. Hybrid PG1159 star Abell 43. 

3.1. Single night observations 

We initially computed Discrete Fourier Transforms (DFTs) of 
the three long light curves acquired using the NOT. We show 
only the DFT from one of the longest runs in Fig.|2] as the other 
DFTs look similar We list the following observed features evi- 
dent from the DFTs of our NOT data: 

- Each pulsation spectrum shows three dominant components, 
which are relatively stable in frequency, but exhibit variable 
amplitudes. 
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Table 2. Identified periodic signals in Abell 43 in the three long 



NOT runs 


Date 


Res. 


Fl 


PI 


Al 


F2 


P2 


A2 


F3 


P3 


A3 


F4 


P4 


A4 


F3/F1 


2005 


(yuHz) 


(y"Hz) 


(s) 


(mma^ 


0"Hz) 


(s) 


(mma) 


(y"Hz) 


(s) 


(mma) 


0"Hz) 


(s) 


(mma) 




06-29 


63 


183.3 


5456 


2.0 








363.9 


2748 


2.3 


416.3 


2402 


0.8 


1.98 


<T 




1.4 


41 


0.1 








1.2 


10 


0.1 


3.6 


21 


0.1 




07-01 


36 


175.5 


5698 


3.1 


321 


3118 


1.0 


366.4 


2729 


2.7 


408.3 


2449 


1.4 


2.08 


cr 




1.2 


38 


0.1 


3.4 


34 


0.1 


0.9 


10 


0.1 


3.4 


15 


0.1 




07-02 


38 


185.2 


5399 


2.1 


325 


3075 


0.8 


368.8 


2711 


1.9 


418.9 


2387 


2.0 


1.99 


cr 




1.6 


47 


0.1 


4.4 


41 


0.1 


1.8 


1.4 


0.1 


3.3 


19 


0.1 





relative amplitude x 10 



- All pulsation spectra have a strong high frequency cut-off 
above 0.5 mHz, beyond which we do not detect any signals. 

- The light curves are non-stationarjO and we have to be cau- 
tious in drawing conclusions using standard tools for spectral 
analysis. 

The dramatic cut off implies that we only find pulsations 
with periods longer than ~2000 s, as predicted both by Quirion 
et al. (I2005I I and Corsico et al. ( 2006) . We determined the fre- 
quencies, amplitudes, and phases of the significant peaks in the 
pulsation spectra using a simultaneous non-linear least squares 
fit {Period04; Lenz & Breger lT005l) . To determine which peaks 
in the DFT constituted real power, we used the method of 
prewhitening in conjunction with the non-linear least squares 
analysis. This involves fitting the dominant component and sub- 
tracting it from the light curve. If we now compute a DFT of 
this resultant light curve, the peak corresponding to the dominant 
mode will be absent. We repeat the procedure for the next high- 
est peak in the pre-whitened DFT, fitting both high-amplitude 
components simultaneously and then subtracting them from the 
light curve. The resultant DFT will not show power due to both 
the highest and second-highest peaks. We continue this process 
as long as there are well-resolved peaks in the DFT. When the 
prewhitened DFT does not show any significant power, then we 
can be sure that we have determined the entire pulsation spec- 
trum using the simultaneous non-linear least squares fit. Fig. [3] 
serves as a demonstration of the prewhitening method. 

We show the results of our analysis for the three longest 
NOT runs in Table |2] that list the significant periods detected 
with a False Alarm ProbabiHty better than 1/100 (Scargle 1982). 
The uncertainties listed in the table are 1 cr errors from the non- 
linear least squares fit. Fig. [T] also includes the synthetic light 
curves constructed from the frequencies, amplitudes and phases 
determined as above. 

It is evident from Table |2] that Fl varies most in frequency 
and amplitude from night to night. These variations may be 
caused by unresolved power We also notice that the ratio of 
the frequencies F3:F1 is very close to two on the first and third 
nights, insinuating a harmonic ratio. The ratio is close to 2.08 on 
the second night and may perhaps be due to Fl being an unre- 
solved multiplet. 

If we compare the periods we detect to those observed ear- 
lier, we find that the period 2473 s discovered by Ciardullo & 
Bond( fT996l l is close to P4, while the period of 5500 s found by 
Schuh et al. ( |2000i) is close to PI. The tw o periods 2600 and 
3035 s, discovered by Vauclair et al. (I2005I I are not too far from 
P3 and P2. The discrepancies may either be due to insufficient 
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' When signals detected in a light curve are not present for its entire 
duration, it is said to be non-stationary in nature. 



300 
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Fig. 3. We show the DFT computed from the combined light 
curve of five runs used in the analysis. The top curve shows the 
window function. The second curve shows the DFT of the light 
curve. The third curve shows the DFT after subtracting out the 
two highest amplitude peaks. The lower four curves show the 
DFT after subtracting three, four, five, and six frequencies re- 
spectively. 

resolution or may perhaps indicate real period changes between 
the different pulsation modes. 

3.2. Combining all useful observations 

To investigate whether the nightly variations in frequencies and 
amplitudes were due to unresolved peaks, we re-analyzed the 
combined light curves with the addition of the APO observations 
from the 10th of July, 2005. Note that we could not utilize the 
APO observations from the 1st of July as bad weather rendered 
them far too noisy for inclusion in our analysis. 

We computed a combined DFT again, and used the 
prewhitening technique in conjunction with the non-linear least 
squares analysis to arrive at the values listed in Table [3] We 
show our final DFT in Fig. [3] along with the different stages 
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Table 3. Identified periodic signals in Abell 43 from the com- 
bined data set 



NGC 7094 NOT 2006-11-29 



Name 


Freq 


Period 


Amplitude 


Noise 




OuHz) 


(s) 


(mma^ 


(mma) 


fl 


164.62 


6075 


0.8 


0.15 


cr 


0.08 


3 


0.1 




f2 


183.76 


5442 


1.9 


0.15 


cr 


0.03 


1.0 


0.1 




f3 


327.70 


3051.6 


0.9 


0.20 


cr 


0.07 


0.7 


0.1 




f4 


363.79 


2748.8 


0.9 


0.20 


cr 


0.08 


0.6 


0.1 




f5 


373.12 


2680.1 


2.0 


0.20 


cr 


0.03 


0.2 


0.1 




f6 


420.23 


2379.7 


1.0 


0.20 


cr 


0.06 


0.4 


0.1 





relative amplitude x 10 



Abell 43 photometry 2005 
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Fig. 4. We show observed light curves, averaged over a sampling 
time of 200 s, compared with synthetic curves constructed from 
parameters in Table [51 



of prewhitening. We find that the power seen earlier as Fl split 
into two individual frequencies f 1 and f2, while the power at F3 
split into f4 and f5. Using the new frequency components, we 
find that f3/f 1 = 1 .99 1 and (f4 -i- f5)/2 f2 = 2.005, giving us two 
ratios of frequencies close to 2. 

Averaging the observations to a sampling time of 200 s, we 
show the observed light curves along with the synthetic light 
curves constructed from our fit in Fig. 21 an inspection of this 
figure reveals that the fit is quite good. This tells us that Abell 43 
is more stable than the pulsating PNNi with shorter periods, in- 
vestigated by Gonzalez Perez et al. (2006). We also notice how 
the light curve can change dramatically on short time scales due 
to beating of the closely spaced modes; this suggests observa- 
tions over a shorter time span compared to the beat cycle may 
give wrong results or show no pulsations at all. 

The shortest period we detect in our data is 2380 s, which 
is quite close to the predicted lower limit of 2600 s for i-l 




69.3 



69.4 

BJD-2454000 



69.5 



Fig. 5. We show the observed light curve of NGC 7094 with a 
simulated light curve based on the periods 34, 50, and 83 min. 



g-mode pulsations induced by the /c-mechanism (Quirion et al. 
120051 and Corsico et al. 120061 1 . The longest period of 6075 s is 
slightly longer than predicted by Qumon et al. (|2005) but within 
the range predicted by Corsico et al. (I2006I I. The non station- 
ary nature of the light curves makes it difficult to draw certain 
conclusions about the frequency ratios close to 2.00. They may 
indicate simple harmonics or non linear resonances as proposed 
by Buchler et al. (fT997l l. 



4. Hybrid PG1159 star NGC 7094: Almost a 
spectroscopic twin of Abell 43 

The hybrid PG 11 59 star NGC 7094 shows almost the same spec- 
troscopic parameters as Abell 43, except that 1% He in Abell 43 
has to be replaced by 1% O (Miksa et al. |2002| l. Quirion et al. 
(120051) propounded that NGC 7094 should pulsate with periods 
in the range of 2550-5413 s. Previous searches by Ciardullo et 
al. ( 11996) and Gonzalez Perez et al. ( |2006l l did not reveal any 
pulsations in the star The latter observation lasted only 5220 s 
and was probably too short to detect pulsations in the predicted 
range. 

In order to scrutinize the variability of NGC 7094, we asked 
for a service run on the NOT, which can be a maximum of 4hr 
long. We obtained a light curve of length 13692 s with ALFOSC 
on the 29th of November, 2006. We used the B filter and a win- 
dow of size Ikx Ik pixels. We chose an exposure time of 30 s 
and the readout time was about the same, giving us a time reso- 
lution of 60 s. Fig. IHshows the light curve we obtained from our 
observations. With the false alarm probability criterion of 1/100 
or better, we identified three peaks with periods between 2000 s 
and 5000 s. A synthetic light curve constructed from these peaks 
is also shown in Fig. |5l Longer observing runs are needed to 
characterize the pulsations better, but this paper serves to estab- 
lish NGC 7094 as a low-amplitude pulsator. 

The range of pulsation periods is shifted toward shorter peri- 
ods relative to Abell 43. This is predicted by Quirion et al. (120051 ) 
and is a consequence of adding a trace of oxygen in their mod- 
els. The observed lower period limit of 2000 s is shorter than the 
theoretical hmit of 2550 s predicted by the same authors. 



5. Conclusions 

For the hybrid PG 1 159 star Abell 43, we find six pulsation pe- 
riods with a range close to that predicted for ^ = 1 g-mode pul- 
sations driven by the /c-mechanism due to partially ionized car- 
bon (Quirion et al. 2005; Corsico et al. 2006). The light curves 
are non stationary, but the pulsation frequencies and amplitudes 
are relatively stable during the 1 1 day span of the observations. 



J.-E. Solheim et al.: Abell43: Longest Period Planetary Nebula Nucleus Variable 



5 



This is different from pulsating PNNi, which showed changes in 
ampUtudes from night to night, and also during the same night 
(Gonzalez Perez et al. 2006). 

We have also detected pulsations in the hybrid PG 1159 
star NGC 7094 which is an approximate spectroscopic twin of 
Abell43. The observed pulsations lie in the range of 2000- 
5000 s, in a band of periods shifted to shorter values than for 
Abell43. This agrees with theory and occurs when a small 
amount of oxygen is added to the model pulsator (Quirion et al. 
12005 ). although the lowest period observed is somewhat shorter 
than predicted for NGC 7094. 
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